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1/12°  Global  HYGHDM :  %  A 
Initial  Development  and  Evaluation 

E.J.  Metzger,  A.J.  Wallcraft,  H.E.  Hurlburt, 
E.P.  Chassignet  and  W.J.  Schmitz,  Jr. 

Ocean  Sciences  Meeting 
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HYCOM  Long-term  Goals  for 
Operational  Ocean  Prediction 

■  1/1 2°  fully  global  ocean  prediction  system 
transitioned  to  NAVO  in  2007 

■  Include  shallow  water,  minimum  depth  5  m 

■  Coupled  sea-ice  model  (LANL  CICE) 

■  Increase  to  1/25°  resolution  globally  by  the 
end  of  the  decade 

■  Optimal  resolution  for  basin-scale 

■  Boundary  conditions  for  coastal  models 
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Global  HYCOM  Configuration 

Horizontal  grid:  1/1 2 “equatorial  resolution 

■  4500  x  3298  grid  points,  6.5  km  spacing  on  average,  3.5  km  at  pole 

Mercator  79°S  to  47°N,  then  Arctic  dipole  patch 

Vertical  coordinate  surfaces:  26-28  for  o0,  32  for  ct2* 

KPP  and  GISS  mixed  layer  models 

Thermodynamic  (energy  loan)  sea-ice  model 

Surface  forcing:  wind  stress,  wind  speed,  thermal  forcing, 
precipitation,  relaxation  to  climatological  SSS 

Monthly  river  runoff  (986  rivers) 

Initialize  from  January  climatology  (GDEM3)  T  and  S,  then  SSS 
relaxation  from  PHC  3.0 

■  No  subsurface  relaxation  to  climatology 
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1/1 2°  Global  HYCOM  snapshot:  SSH  and  ice  (gray) 


■  Running  at  NAVO  under  DoD  Challenge 

■  190K  CPU  hrs/model  year  on  784  CPUs 


6  Dec 2005  ■  7  2  TB/model  year  for  daily  3-D  output 


1/1 2°  Global  HYCOM  Experiments 

■  ECMWF  Reanalysis  (ERA15)  climatological 
wind  and  thermal  forcing 

■  Annual  bias  corrections  to  air  temperature 
(ERA40),  radiative  fluxes  (ISCCP)  and 
precipitation  (GPCP) 

■  o0  simulations: 

■  26-layers  with  KPP  for  6  model  years 

■  28-layers  with  GISS  for  3  model  years 

■  o2*  simulations: 

■  32-layers  with  GISS  for  9  model  years 
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Initial  1/1 2°  Global  HYCOM  a0  Simulation 

■  Major  shortcomings: 

■  Poor  simulation  of  both  Gulf  Stream  and  Kuroshio 

■  Poor  representation  of  tropical  current  systems 

■  Unrealistic  transport  at  key  locations: 

■  Florida  Straits  (23  Sv  vs.  -32  Sv)  [simulated  vs.  observed] 

■  Drake  Passage  (91  Sv  vs.  -134  Sv) 

■  Pacific  to  Indian  Ocean  Throughflow  (22  Sv  vs.  -10  Sv) 


Improved  1/1 2°  Global  HYCOM  ct0  Simulation 

■  Modifications 

■  Added  two  layers  (26  ->  28)  and  changed  layer  structure 

■  Increased  eddy  viscosity:  A  =  30  m2/s  constant  everywhere 

■  Increased  Smagorinsky  diffusion:  .05  ->  .1 

■  KPP  GISS 
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Improvements  in  Tropical  Pacific  Current  Structure 


Observations 


Original  26-layer  o0  Modified  28-layer  o0 
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1/1 2 “Global  HYCOM  Mean  Gulf 


Stream  And  Kuroshio  Pathways 


1/1 2°  Global  HYCOM  o2*  Simulation 

The  Atlantic  subpolar  gyre  generally  looks  good 
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Global  SSH  Variability 
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Oct  92  -  Nov  98  SSH 
variability  based  on  T/P, 
ERS-1  and  ERS-2  altimeters 
(from  Collecte,  Localisation, 
Satellites  (CLS)) 


SSH  variability  from  1/12° 
global  HYCOM  a2*  with 
climatological  wind  and 
thermal  forcing 
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Transport  (Sv) 


Atlantic  Meridional  Overturning  Circulation 

a0  vs.  a2* 
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DEPTH  (m) 


Atlantic  Meridional  Overturning  Circulation 
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Velocity  Cross-section  Along  Luzon  Strait 

Sb-ADCP  data  (top)  vs.  1/1 2°  global  HYCOM  (bottom)  in  the  upper  300  m 
Section  along  21  °N  between  1 18.5°E  and  124.0°E 


Sb-ADCP  data  from  Liang  et  al.  (2003,  DSR  Pt.  II) 

Mean  from  HYCOM  with  ERA15  wind  and  thermal  forcing 

No  ocean  data  assimilation  in  HYCOM  16 


Vertical  Structure  in  the  Equatorial  Pacific 


Observations 
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:  Zonal  temperature  along  the  equator 
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Deep  Flow  Through  Samoan  Passage 


Observed  mean  northward 
transport  below  4000  m  =  6.0  Sv 
Rudnick  (1997,JGR) 


HYCOM  ct2*  in  layers  27-32  =  9.3  Sv 
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Transport  Comparisons  at  Key  Locations 


Section 

Obs. 

Orig.  a0 

Mod.  o0 

Orig.  a2* 

Mod.  a2* 

PCM-1 

23 
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25.5 
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Denmark  Strait 
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Drake  Passage 

134 

91.3 

96.2 

152.4 

146.4 

Modified  o2*  experiment  uses  a  new  topography  with  sill  depth  refinements 
in  the  IAS,  Indonesian  Seas,  etc. 
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Future  Work  (FY06) 

■Ten  May  2001  -  June  2002  assimilative  runs  in  FY06 

■  Time  period  with  three  satellite  altimeters 

■  Five  with  bi-weekly  30-day  forecasts 

■  Near  real-time  nowcast/forecast  starting  in  mid-FY06 

■  Interannual  non-assimilative  case: 

■  1995-present  using  NOGAPS 

■  Coupling  with  LANL  CICE  via  ESMF 
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